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SUMMARY 

Samples of four silica formers: single crystal SIC, sintered a-S1C, 

reaction sintered SI 3 N 4 and polycrystalllne M 0 SI 2 , were subjected to a Mach 1 
jet fuel burner for 1 hr, at a sample temperature of 1375 °C (2500 °F). Two 
phenomena were Identified which may be deleterious to a gas turbine appli- 
cation of these materials. The glass layer formed on the M 0 SI 2 deformed 
appreciably under the aerodynamic load. A scale developed on the samples of 
the other materials which consisted of particulate matter from the gas stream 
entrapped In a S10 2 matrix. 


INTRODUCTION 

Silicon based materials are under consideration for a number of applica- 
tions in aircraft and automotive gas turbine engines. The majority of these 
are monolithic components which benefit from the excellent high temperature 
strength of some of these materials, eg. sintered alpha silicon carbide and 
reaction bonded silicon nitride. Another possible application Is employing 
these materials In turbine tip gas path seals. The structural requirement of 
this application Is minimal, however dimensional stability Is crucial. 

The silicon based material are unstable In the presence of oxygen at the 
high temperatures experienced In the engine. Protection from gross reaction 
Is afforded when an adherent film of SIO 2 forms on the surface during oxida- 
tion. This film serves as a barrier across which mass transport must occur 
for further oxidation, leading to parabolic scale growth kinetics. The class 
of materials which displays this behavior Is termed "silica formers." 

It has been postulated that one-dimensional growth of an oxide film with 
a positive volume change, and parabolic kinetics on the surface of a turbine 
tip seal can be favorable (ref. 1). An additional requirement for this Is 
that the oxide film Is not deformable under the frictional forces produced by 
gas flow through the engine. 

The materials studied were; molybdenum dlslllclde, single crystal silicon 
carbide, sintered alpha silicon carbide, and reaction bonded silicon nitride. 
It Is well documented (e.g. ref. 2) that under a static air environment, at 
temperatures In excess of 1500 °C (2730 °F), these materials oxidize to pro- 
duce an adherent scale of vitreous silica. The purpose of this work was to 
examine the morphology of scales produced when oxidation took place In the 
high velocity effluent of a (Mach 1) jet fuel burner. 





EXPERIMENTAL 


The sample sources, geometries, and surface preparations are detailed In 
table I. 

The burner rig used was the NASA Lewis High Temperature Erosion Rig which 
Is described In detail In reference 3. The essential components of the rig 
are a jet fuel combustor, a Mach 1 nozzle, and a sample holder, see figure 1. 
The rig also has provision for Introducing grit particles to the gas stream, 
however this was not used. 

The samples were placed In sample holder 50 mm (2 Inch) from the nozzle 
exit slightly Inclined, 15°, from perpendicular. Temperatures were measured 
using pyrometers. The same rig parameters (fuel pressure, air pressure, etc.) 
were used In all cases. Some variation was found In the temperature measure- 
ments with readings between 1325 and 1390 °C (2420 and 2535 °F). The diameter 
of the flame on the sample was 8 mm, (0.3 Inch) outside of this a steep thermal 
gradient existed. Sighting of the pyrometer may have contributed to the scat- 
ter. The time of exposure to the flame was 1 hour In all cases. 


RESULTS 

Molybdenum Dlslllclde 

Two types of surface preparation were used on the M 0 SI 2 rods: as ground 

and static preoxidation at 1300 °C (2340 °F) for 24 hr. Both types were run 
In the Mach 1 burner rig for 1 hr at 1375 °C (2510 °F). The post test condi- 
tion of the rods Is shown In figures 2 and 3. From these photographs It Is 
apparent that the SIO 2 produced through oxidation was physically blown from 
the leading edge and accumulated on the leeward side of the rod. This leaves 
the MoS 1 2 unprotected and allows the reaction to proceed as evidenced by the 
extensive pitting on the rod surface. The sample given the static oxidation 
appears to have less severe corrosion than the as-ground sample, however It Is 
clear from the photographs that the preoxidation did not affect the character 
of the reaction. 


Single Crystal Silicon Carbide 

The single crystal silicon carbide was exposed to the burner In an as- 
polished state. The behavior was qualitatively different from that of the 
M 0 SI 2 . Scaling, or a buildup of material on the SIC surface occurred directly 
under the location of the flame, see figure 4, while the rest of the surface 
displayed Interference colors typical of passive oxidation. A typical cross 
section of the scale material Is shown In figure 5. The thickness (50-100 ym) 
Is far In excess of what Is expected from the formation of SIO 2 during pas- 
sive oxidation («0.2 ym) for 1 hr at 1375 °C (2510 °F). 

Two sources of unanticipated sources of particulates were Identified after 
the experiments were run. The first was residual alumina grit particles from 
previous operation of the rig In the erosion mode. Evidently a certain frac- 
tion of the grit particles which had accumulated In the plumbing of the rig 
became entrained In the gas stream during our experiments. It should be noted 
that no erosion due to particle Impact has been observed for either oxide or 
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metal specimens under these conditions. This Implies that the number density 
of these residual particles Is small. 

The second source of particulates was the sample holder Itself. As Illu- 
strated In figure 11, the stainless steel clamping block used to hold the spe- 
cimen had a corner which protruded very close to the main gas stream. Some 
erosion of metal from this corner was observed. These erosion products would 
be expected to become entrained In the gas stream. 

The microstructure of the scale Is complex. There Is an Intermlttant 
contact between the scale and the carbide. The scale also appears cracked and 
contains Inclusions. Energy dispersive x-ray analysis revealed the presence 
of SI, A1 , Fe, and T1. Using an electron microprobe, element mapping was done 
near the carbide/scale Interface. Haps for SI, Al, Fe, and 0 are shown In 
figures 4 to 7 respectively. It Is clear from these that silicon and aluminum 
are present as oxides while Iron Is present as a metallic species. A second, 
and striking, feature Is that the silica phase Is the minority phase by volume 
but It forms a continuous matrix which contains Inclusions of alumina and Iron. 
Reference to the micrographs also Indicates that the scale thickness Is several 
times the diameter of the Inclusions. These micrographs therefore suggest a 
mechanism for scale growth by the adhesion of particles In a "glue" of silica. 
This can give rise to the formation of scales which are many times thicker 
than what Is predicted on the basis of oxidation alone. 


Sintered Alpha Silicon Carbide 

Samples of both as-fabricated and as-fabricated plus static oxidation at 
1300 °C (2370 °F) for 24 hr were exposed to the burner. The preoxidation had 
no effect on the observed result. The scale growth on the sintered alpha 
silicon carbide appeared qualitatively the same as the single crystal silicon 
carbide except that the boundary of the scaled region Is less distinct, see 
figure 12. Attempts at a transverse section were unsuccessful. 

It was also observed, under high magnification, that the scale contained 
a high density of bubbles, see figure 13, which Is typical for this material 
(ref. 4). 


Reaction Bonded Silicon Nitride 

One sample of silicon nitride In as-fabricated form was exposed to the 
burner. A scale formed on the block which visually appeared to be similar to 
that formed on the silicon carbide. It was not Investigated on the micro- 
scopic level. 


DISCUSSION 

Identification has been made of two phenomena which may affect the suit- 
ability of the use of silica formers In gas turbines. The first was the defor- 
mation of the oxide scale film under the aerodynamic load. This Is highly 
undesirable because the existence of a continuous mass transport barrier Is 
fundamental to passive behavior. It Is notable that the S10 2 film on the MoS1 2 
was observed to be much more deformable, Inferentlally of a lower viscosity. 
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than that on the silicon carbide. The reason for this Is unclear, It Is 
possible to speculate that the presence of molybdenum reducing the viscosity 
of the scale. 

In the appendix a simple model Is developed to determine the deformation 
and mass flow of a viscous glass film due to aerodynamic friction over an 
oxidized rod. The model was applied using molybdenum dlslllclde parabolic 
rate constants (ref. 5) at various stagnation point (on the rod) oxide layer 
thicknesses and viscosities. The model results are shown In figure 14. The 
observed mass flow suggests a viscosity several orders of magnitude lower than 
that reported for vitreous silica at this temperature (viscosity 
*1.0xl0 9 N«s/m 2 at 1450 °C). 

The severity of this problem might also be enhanced during long time 
exposures when Impurities In the fuels (such as Ca, Cu, Fe, and Pb) may 
deposit on and reduce the viscosity of the scale. 

The second phenomenon observed was the adhesion of debris within a matrix 
of SIO 2 somewhat analogus to bricks and mortar. The SIO 2 was not the major 
phase by volume but served the role as the binder. This would be a problem In 
a seal application If debris collects on the seal surface producing high spots 
that lead to rubs against the turbine blades. It would also be undesirable to 
have the shape and surface roughness of airfoils made from silica formers modi- 
fied by the deposition and adherence of debris. While the types of debris 
observed were particular to the test set up It Is anticipated that any type of 
particulate Ingested or produced In an engine would likely adhere. 


CONCLUSIONS 

These experiments were somewhat preliminary In nature and further experi- 
mentation under more controlled circumstances Is desired, however, It may be 
concluded that: 

1. The SIO 2 film formed on the M 0 SI 2 rods deformed appreciably under 

the aerodynamic load of the high velocity effluent of a Mach 1 jet fuel burner. 

2. The scale formed on single crystal silicon carbide was composed of 
metal and oxide particulates embedded In matrix of SIO 2 . These scales were 
several particle diameters thick and many times thicker than expected from 
oxidation alone. 

3. The scales formed on sintered alpha silicon carbide and reaction bonded 
silicon nitrides appeared similar to the scale formed on the single crystal 
silicon carbide. 
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APPENDIX: EFFECT OF HIGH TEMPERATURE COMBUSTION FLOW OVER 

AN OXIDIZED ROD 

A simple flow model is developed to determine the deformation of a viscous 
glass film due to aerodynamic friction at high temperatures as a function of 
position on the periphery of the oxidized rod. The overall model Is shown In 
figure (15) with the parameters of Interest Indicated. 

The assumptions used for this analysis are as follows: 

(1) constant gas temperature equal to 1700 K. 

(2) Constant gas velocity equal to sonic velocity for air at 1700 K. 

(3) Linear pressure drop (constant pressure gradient) around the peri- 
phery from the stagnation point to the separation point (ref. 6). Pressure 
was assumed to be the same In the highly viscous layer as In the gas. 

(4) Shear stress due to the hot gas acting on the oxide layer was constant 
and was found from a potential flow solution to the Navler-Stokes equations 
for cylindrical geometry (ref. 7). 

(5) The rate constant for oxide layer growth was assumed to be constant and 
given In reference 8. 


Analysis 

First the continuity equation Is applied to the model and Is shown In 
figure 16. The change In mass flow, Am, is given by the following equation: 


Am = (RA0) (AZ) pM 

where 

R radius of the rod 

A0 angular Increment 

AZ depth (length) of rod 

p density of viscous layer (oxide) 

M rate of flow of rod material to oxide layer 

The mass flow rate, M, of the rod to oxide layer Is given by 

K K 

M = JB = 
i .2 

6p 


(Al) 


( A2) 


K p parabolic rate constant from reference A3 
6 viscous layer thickness 
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The rate constant K p was taken from reference 8 and Its value used In 
this analysis was 

K . 4.8x1 0 - ^ 

P m 4 . S 

which was the experimental value found In reference 8 for molybdenum 
dlslllclde at 1645 K and partial pressure of oxygen of 14 torr. 

Substituting equation (A2) Into (Al) gives: 

(iz) (iS) ■ ** (?) < A 

letting Am be the change In mass flow rate unit depth and letting A0 
approach zero gives: 



( A4) 


This expression determines the rate of flow of rod material with respect to 
angular position and Is dependent on the oxide layer thickness. 


Next the momentum equation Is needed to determine the loading effect on 
the oxide layer thickness. Figure 17 shows the model used for this 
formulation. 


The forces that are on the oxide layer are shown In figure 17. For a 
nondeformable control volume In steady flow the sum of the forces Is given by 
(ref. 6): 


Z F ■ X.s. pV(V ' ft) " <fl5 > 

For low Reynolds number flow the flux terms on the right hand side of 
equation (A5) can be neglected. From the diagram this gives 

- (AP6L) + t 0 ( R A6) L = tRL A 0 (A6) 


or: 


- (AP5) + Tq(R A0) = tR A0 


Dividing through by A0 



+ T R 
o 


tR 


where 


AP 

A0 


( A7) 


(A8) 


pressure gradient 
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t 0 ,t shear stress gas to oxide, oxide to rod 

The gas pressure was assumed to decrease linearly as a function of 0, 
from the stagnation point to the separation point (fig. 15). The pressure was 
assumed to behave by Bernoulli's Theorem given by reference 6 as: 

P 1 + \ P V 2 = C 1 (A9) 


which means the total pressure Is constant or 

W - C 1 - C 2 e i 


ae. 


= - c. 


(A9a) 


where P^ , pressure and velocity at position 0^. 

The constants, C-j and C 2 were found by using the property data of Jet A 
fuel burned In air found In reference 9, and the assumed sonic velocity of the 
high temperature gas. For this analysis C-j and C 2 were found from 

P = 0.2122 Kg/m 3 , 

V = 826 m/s, and 

P = 1.0133xl0 5 N/m 2 ; 

C ] = 1.737x10 s N/m 2 

C 2 = 3.7878X10 4 N/(m 2 *rad1an) 

The next parameter needed In equation (A8) was, t 0 , or the shear 
stress acting on the gas-oxide surface. An Integrated average value was found 
from figure 18. This figure plots the nondlmenslonal shear stress versus 
nondlmenslonal position for various types of bodies. These curves were 
developed from a potential flow solution of the Navler-Stokes equations 
(ref. 7) for various shaped bodies. Using the parabolic shape of the circular 
cylinder plot^an Integrated average nondlmenslonal shear stress was found. 

This value, x 0 , was found to be equal to 4.066. This value along with the 
gas properties provided the average shear stress x 0 . 



U free stream gas velocity; (U = V ) 
00 3 J ’ ' <x> sep 

p gas density 
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V 


kinematic viscosity of the gas 


d diameter of the rod 

Solving for t 0 with the appropriate values gives: t 0 = 2106.7 N/m? per 

unit depth. 

To evaluate- the shear stress between the viscous oxide layer and the rod, 
t, It Is necessary to evaluate the flow of the oxide layer first. 

= pijl V e (All) 

m^ mass flow rate of oxide layer at 1 th position 

p density of viscous layer 

5^ oxide layer thickness at position 1 

L length of rod exposed to gas flow 

Ve average velocity across oxide layer thickness 

If It Is assumed that the velocity profile In the oxide Is parabolic, then the 
velocity as a function of position In the oxide Is given by: 

V e (y) = A + By + Cy 2 (A12) 

where y Is the coordinate In the normal direction from the rod surface. The 
following boundary conditions are Imposed: 

Boundary Condition 1: V 0 = 0 at y = 0 (at the rod-oxide Interface) 

aV e 

Boundary Condition 2: = 0 at y = & 

Using condition (1) In equation (A12) gives 

A = 0 

using condition (2) In equation (A12) gives: 


B = -2 C5 


Substituting A and B Into equation (A12) 

V Q (y) = -2C6y + Cy 2 (A13) 

For the continuity equation the average velocity, V 0 , Is needed and this Is 
given by reference 6 to be: 


v e (y) 



(A14) 
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Substituting equation (A13) Into (A14) and evaluating gives the following 
expression for average velocity: 



(A15) 


Inserting equation (A15) Into equation (All) and solving for C gives: 


C 


3m 

2 P & 3 


(A16) 


Now using equation (A13) and the relationship between shear stress and 
velocity profile: 


T 



y=o 


(shear stress at the Interface) 


(A17) 


rf = -2C6 (A18) 

dy I y=0 

Substituting equation (A18) Into (A17) the shear stress at the oxide-rod 
Interface Is: 



(A19) 


v oxide layer viscosity 
m flow rate of oxide per unit length 
p oxide layer density 

6 oxide layer thickness 


Now the shear stress between the rod and the oxide layer Is known In 
terms of the oxide layer flow rate and the oxide layer thickness. Now 
equations (A19), (A10), and (A9a) can be Inserted Into equation (A8) and 
solved for "i", the oxide layer thickness. This equation then becomes: 


& 
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(A20) 


Numerical Method 

The numerical model used Is shown In figure 19. The procedure used Is as 
follows: 

(1) Assume an Initial layer thickness 6 0 at the first step, set m 0 = 0. 
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(2) Use the Initial guess 6 0 with equation (A4) and a fourth-order 
Runge-Kutta approximation to find the mass flow rate per unit depth 

at the 1™ position. 

(3) Using Si just found, substitute Into equation (A20). Note: 


(4) 

(5) 

( 6 ) 


m 


1+1 


ii^ + diti/de 


61 


Solving the third degree equation by Newton^s Method for 6^ + -|. 
This Is now used at the next position for m In equation (A4). 
Repeat the process until the separation point Is reached. 


Numerical Results 

The numerical results of the above analysis are shown on figure 14 as 
presented In the discussion. This figure Is a summary of varying the Initial 
layer thickness 6 0 and the oxide layer viscosity on the mass flow rate per 
unit depth at the separation point. At high viscosity the behavior of the 
mass flow Is dominated by the oxide thickness, 6 0 , at the stagnation point. 
However, as the oxide layer viscosity decreases the effect of a thin Initial 
layer Is minimal on the mass flow rate at the separation point on the rod. 
This Is evident as the plotted mass flow rate per unit depth Is Invariant for 
two orders of magnitude difference In the Initial oxide layer, 6 0 , thickness. 
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TABLE I. - SPECIMEN TYPE, SOURCE, AND GEOMETRY 


Sample I.D. 

Source 

Burner rig 
sample geometry 

a Molybdenum dlslllclde 
(2 specimen) 

b $1ngle crystal silicon carbide 

Wintered alpha silicon carbide 
(2 specimen) 

Reaction bonded silicon nitride 

Kanthal Corp. 

Carborundum Co. , 

Acheson furnace 0 
Carborundum Co. c 

Air Research Casting Co. 

rods, 5 mm (0.2 Inch) 
diameter 

crystals, roughly 10x10x3 mm 
(0.4x0.4x0.12 Inch) 
quarters of a tube, 

18 mm (0.7 Inch) o.d. 
blocks, 13x25x50 mm 
(0. 5x1x2 Inch) 


a $urface preparation: First specimen etched with HF (25%) for 45 min, finished ground 

with 325 mech SIC. Second specimen same as first then statlcly oxidized In air at 1300 °C 
(2370 °F) for 24 hr. 

^Surface preparation polished to scratch free surface with successive grades of diamond 
paste down to 1/4 ym. 

c Suppl1ed by J.W. Halloran Case Western Reserve University. 

d Surface preparation: first specimen as fabricated, second specimen statlcly oxidized In 
air at 1300 °C (2370 *F) for 24 hr. 

e Surface preparation; as fabricated. 
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Figure 5. - SEM micrograph of the cross section of the builtup material using backscattered electrons. 







Figure 9. - Iron elemental distribution around carbide/scale interface. 




Figure 10. - Oxygen elemental distribution around carbide/scale interface. 



Figure 11. - Specimen holding fixture showing 
the orientation of test hardware in the jet 
fuel combustion products. 
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Marne direction 


Figure 12. - SEM micrograph of the scale built up on the sintered alpha silicon carbide. 




Figure 14* ~ Effect of initial layer thickness on mass flow 
per unit depth at the seperation point for a wide range 
of viscosity of oxide layer* {M 0 Si 2 rod). 



Figure 15. - Model of system analyzed. 


r Viscous 
I oxide layer 



Figure 16. - Continuity equation model. 




Non dimensional position, x/L 


Figure 18. - Shear stress as a function of position for various geometry. From Fluid 
Dynamics Handbook by Streeter. For circular cylinder: X - 0| ; L = 9 m at X/L - 
. 609 separation point a/b - aspect ratio for elliptic bodies. 



Figure 19. - Numerical model. 
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